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(N 

(**) I' Abstract. We introduce n-parameter M d -valued Brownian-time Brownian 

^^ sheet (BTBS): a Brownian sheet where each "time" parameter is replaced with 

the modulus of an independent Brownian motion. We then connect BTBS to 
'\^J , a new system of n linear, fourth order, and interacting PDEs and to a corre- 

Cn ■ sponding fourth order interacting nonlinear PDE. The coupling phenomenon 

is a result of the interaction between the Brownian sheet, through its variance, 
and the Brownian motions in the BTBS; and it leads to an intricate, intrigu- 
Q—( ing, and random field generalization of our earlier Brownian-time-processes 

P^ ■ (BTPs) connection to fourth order linear PDEs. Our BTBS does not belong 

to the classical theory of random fields; and to prove our new PDEs con- 
nections, we generalize our BTP approach in [5] [4] and we mix it with the 
Brownian sheet connection to a linear PDE system, which we also give along 
with its corresponding nonlinear second order PDE and a 2n-th order linear 
PDE that we also connect to Brownian sheet. In addition, we introduce the 
n-parameter d-dimensional linear Kuramoto-Sivashinsky (KS) sheet kernel (or 
#^1 , "transition density"); and we link it to an intimately connected system of 

^. i new linear Kuramoto-Sivashinsky-variant interacting PDEs, generalizing our 

*vj ' earlier one parameter imaginary-Brownian-time-Brownian-angle kernel and its 

^\ ' connection to a linear KS PDE. The interactions here mean that our PDEs 

r^\ , systems are to be solved for a family of functions, a feature shared with some 

jyv ■ well known fluids dynamics models. The interacting PDEs connections estab- 

• ' lished here open up another new fundamental front in the rapidly growing field 

of iterated-type processes and their connections to both new and important 

higher order PDEs and to some equivalent fractional Cauchy problems. We 

t_^ , connect the BTBS fourth order interacting PDEs system given here with an 

interacting fractional PDE system and further study it in another article. 
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1. Introduction and results statements 

1.1. A brief orientation and motivation. In 0|4], we introduced a large class 
of interesting iterated processes that we called Brownian-time processes (BTPs Jj — 
stochastic processes in which the regular clock t is replaced with the Brownian 
clock |-B(£)|, where B is a one-dimensional Brownian motion starting at — and 
we connected them to novel fourth order PDEs. The simplest such BTP is the 
Brownian-time Brownian motion (BTBM) X| := {X x (\B(t)\);0 < t < oo} where 
X x is an independent Brownian motion starting at x £ Mr; and we showed in [5] 
that the associated PDE, for X^ as well as for a large class of related BTPs, is 

(ii) l§i u{t > x) = ^T + l A * u{t ' x); t>0 > x€Rd > 

\u(0,x) = f(x)\ xeM d . 

I.e., if u(t,x) — E/(X|g(i)), then u is the solution to the PDE (|1.1[) for an appro- 
priate class of initial functions / (see also Allouba [I], DeBlassie [IB], and Nane 
[2"3"]). The novelty in the PDE (|l.l|l is the A x f(x)/y/8 / Kt term which expresses the 
BTP memory-preserving (non-Markovian) property. Such PDEs connections are 
fascinating from many perspectives, here are some: 

(1) probabilistically (|1.1|) is the PDE solved by "running" the intrinsically- 
interesting BTP, giving a dynamic description of the average of the modu- 
lated process f(X%); 

(2) the term A x f(x)/\/8irt results in an interesting unconventional interaction 
between the initial data / and the solution u (this feature is under investi- 
gation for its potential applications to financial models); 

(3) along with Zheng, we also initiated in [5] the link between Brownian-time 
processes and fractional time derivatives by way of the BTP half-derivative 



A key of frequently used acronyms is provided in Appendix [B] for the reader's convenience. 
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generator — an implicit equivalence of (|1.1|) to temporally-fractional PDEs 
that was later made explicit and extended in several directions in the nice 
works of Nane, Meerschaert, Baeumer, Vellaisamy, Orsingher, and Beghin 

[si 122 naming; and 

(4) in [3] , we adapted our BTPs construction in [H [5] and their PDEs connec- 
tions and gave an explicit representation of the solution to a linearization of 
the celebrated Kuramoto-Sivashinsky (KS) PDE in all spatial dimensions. 
Besides giving a probabilistic Brownian-time-motivated and explicit repre- 
sentation of this important equation, this achievement is noteworthy since 
the existence of the KS semigroup is still unsettled for dimensions larger 
than 1 via traditional analytical methods, precluding the possibility to have 
definitive results using semigroup methods in higher dimensions. We will 
have more to say about that in upcoming work. 

In this article, we introduce a natural random field generalization of our Brownian- 
time Brownian motion in J5[ 0], which we call Brownian-time Brownian sheet 
(or BTBS); and we also introduce an interesting multi-time-parameter Kuramoto- 
Sivashinsky sheet kernel that generalizes our single-time parameter KS kernel in 
[U [2J Q]. We then link them both to intriguing fourth order interacting PDEs. 
These new PDEs connections add novel and intricate new features not present in 
their one parameter (n = 1) d-dimensional (d > 1) counterparts in [3]. Our 
interest in these fourth order systems of interacting PDEs is at least threefold (1) 
their connections to the intrinsically interesting BTBS random field in Subsection 
11.21 (2) their novel interaction terms which are fascinating, characteristic of the 
multi-parameter Brownian-time setting, and intricate nontrivial generalizations of 
the important one parameter case; and (3) interacting PDEs are important in tur- 
bulence and fluids dynamics models as in the venerable second order Navier-Stokcs 
equations in which the interaction between pressure, velocity, and an external force 
(e.g., gravity) plays a crucial role. In that light, we believe that the interactions 
in the PDEs of this article — which are probabilistically-naturally-derived multipa- 
rameter generalizations of our single-time parameter, fourth order, and turbulent 
pattern formation equations — deserve more investigation. We state and further de- 
scribe the two main results of this paper (Theorem 11.11 and Theorem 11.21) in the 
following two subsections. 

1.2. Brownian-time Brownian sheet and its fourth order interacting PDEs 
system. An interesting and natural random field generalization of our BTBM 
in [5j H] is the Brownian-time Brownian sheet, which we now introduce. Let 
B^ l \ . . . , £?(") be n independent copies of a standard one-dimensional Brownian 
motion starting at and independent of an n-parameter (n > 1) Revalued Brow- 
nian sheet 

W° = {W°(t) = (W?(t), . . . , Wj(t)) ;t = (ti, ... ,tn) G R';} , 

"starting" at <E K d under F—P{W°(t) = 0} = 1 for t E <9R™ (see Notation [HE])— 
all defined on a probability space ($7,J^",P). Of course, the Brownian sheet co- 
ordinates {W° (t); t > 0} , . . . , {W2(t);t > 0} are assumed independent. For any 
x = (xi,...,Xd) G M d , let 

W x :=W° + x = {W x (t) = (WT(t), . . . ,W**(t)) ;t = (t u . . . ,t n ) G R^} . 
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We define the n-parameter Revalued Brownian-time Brownian sheet (BTBS) based 
on W x and B^\ ..., B<- n \ and starting at x £ R d , by 



(1.2) 



Clearly, P 



'bW....j(« 



w- 



B^\...,B( 



(t):=W x {\B^(t 1 )\,...,\B^(t n )\) 
(n , (t) = x] = 1 for t € SRIfB 



t G 



Our first main result (Theorem 11.11 below) gives a novel system of linear fourth 
order interacting PDEs (| 1 -4[) and a corresponding nonlinear fourth order interacting 
PDE (|1.5[) that are solved by running the Brownian-time Brownian sheet random 
field in (|1.2p . When n > 1, this new BTBS-PDE connection describes the dynamics 

of the average of the modulated Brownian-time Brownian sheet / ( W„ (1) B{rl) ) 
as it interacts with the averages of the intimately related random fields given by 

(lLeN„\{j} \ B(l) |) / ( W s(D : ...,B(")) . for ea - ch J G N„ := {1, . . . , n}. Before stating 
our main result, it is helpful to adopt some simplifying notational conventions. 



}; and we denote by 



Notation 1.1. We always use the notation N„ := {1, . . . , n 

o 

RI? and <9R™ the interior and boundary of R™ , respectively. From this point on, we 
will alternate freely between ut, u(t), and u(t±, . . . , t n ) for typesetting convenience 
and ease of exposition. Moreover, ty = (U)ieN n \{j} G R™ -1 will denote the (n — 1)- 
tuple point in R™ _1 that is obtained from t = (ti,... ,£„) £ R™ by removing the 
j'-th variable, tj. The notations ut, u(t), and utJtj) will all mean the n-parameter 
function u : R" —¥ R d evaluated at t = (ti, ... ,t n ), and we use UtJtj) whenever 
we need to highlight the variation in the variable tj for clarity. The same comment 
applies for the notations u{t,x) and UtJtj,x). 

We are now ready to state our first main result. 

Theorem 1.1 (Brownian-time Brownian sheet: fourth order interacting PDEs 
connections). Let {W„ (1) „ ( „)(t);t £ R™} be an n-parameter M. d -valued BTBS 

based on a Brownian sheet W x and Brownian motions {B^>\ . and starting on 
<9R" at x £ R d on (17, J?, P). Let f : R d — > R be bounded and measurable such that 
all second order partial derivatives d^f are bounded and Holder continuous, with 
exponent < a < 1, for 1 < k, I < d. If 



u(t,x)=E P [f(w* B(1) _ B(n) (t) 



(1.3) 



ty( j \t,x)=E w 



n \b {i) (u)\) f (^(i,,... )fl (.) 

^eN„\{j} / 



(t) 



; j e N„, 



for(t,x) £ W£ xR d ; then the family {u(t,x), tyW(t,x), . . . , fy^{t, x)} is a solu- 
tion to both 



We note that in |14l Burdzy discussed a process intimately related to the case n = 2 of our 
BTBS in the context of a question he posed regarding the path properties of such a process. 
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W 



(1.4) 



(ii) 



(1.5) 



(1.6) 



the system of fourth order interacting Linear PDEs: 

AJ(i) + ^A^W; t G K'j ,x G W l , 



du 



lLeN„\{i}** . ., - ' 



2 4 -™7r"i 1 



(&)u(t,s) = f(x); 

(c) Wt»(t,x) = 0; 

(d)^W(t,x) = [a 6N „\ {j} ^/(^); 



t G SKIjUzG R d , 

t y e dK™ _1 ,xe IE 

t, = 0, a; G M d , 



/or j G N„ ; and 

t/ie corresponding fourth order nonlinear interacting PDE: 



TT — 
11 at,- 



"'"^''(Ajwr + in^' 



2 4 



2 7r" 2 



n—n^^-n 



(fci,...,fc„)6S„ieN n 
together with the temporal boundary conditions (6), (c), and (d) m (|1.4p . 



whe 



Tkjjfax) 




2 4 -"7r"t 1 



-A^W; 



A x f(x): i^ljeN, 

fej = 2,j G N n , 



id w/i 



and where 



S„ = i {k x ,...,k n );ki,...,k n e {l,2},n + l<J2 k 3 <2n-l >, n > 2. 



The intriguing interaction (or coupling) in Theorem 11.11 between it and W^', 
j = l,...,n, shows that the PDEs in (I1.4J) and (|1.5[) are nontrivial intricate gen- 
eralizations of our BTP-PDE connections in [51 3] (the case n = 1). This pre- 
cise coupling phenomenon is caused by the interaction between the Brownian-timcs 
\B^ I , . . . , \B( n > | and the variance of the Brownian sheet W x through cjj-Var (W x (t)) 



J| ieN > r ,i t%. The nonlinear feature of the BTBS-PDE connection and its corre- 
sponding linear system connection are inherited from the nonlinear Brownian sheet 
PDE together with its corresponding linear system of PDEs, which we give below in 
Lemma 1 1.1 1 However, in addition to the new coupling above, and due to the itera- 
tive nature of BTBS, those PDEs in (|1.4[) and (|1 .5|) are fourth order versus the much 
simpler second order Brownian sheet nonlinear PDE and linear system of PDEs. On 
the other hand, when n = 1, our BTBS becomes a Brownian-time Brownian motion 
(see [5J 15), u = &&>, and (JOJ and (OJ) reduce to our fourth order BTP-PDE 



|BW( t . 



1 



Tkj.jitjx) = since 



(ED ((0.2) in 0). Note that if n = 1 then rU N „\ {j} U = U iG N n \{j} 

since N„ \ {j} = 4> (the empty set), and T,(k u ...,k n )e§ n Hje 

S n = <fi. Moreover, the interaction in the (|1.4|1 means that the PDEs are to be 

solved for a family of functions \u, ty£i\ as opposed to just one function u as 

in either the Brownian sheet case (|1.7[) or the one-parameter case (|l.l|l . 

As with its one-parameter BTP counterpart, the BTBS is memory preserving 
as is indicated by the inclusion of the Laplacian of the initial function A x f{x) in 
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the BTBS PDEs in (fT4| and (|L5]l . Also, just as BTPs are not classical (not semi- 
martingales, not Markovian, and not Gaussian), BTBS is not a classical random 
field. To prove Theorem 11.11 we generalize and mix our Brownian-time approach 
from OH] with the Brownian-sheet-PDEs connection which we give in Lemma ll.ll 
below. 

Theorem 11.11 gives intriguing connections between the Brownian-time Brownian 
sheet and a novel nonlinear fourth order and interacting PDE (|1.5|) and its corre- 
sponding interacting system (|1.4[) . This connection opens up another new funda- 
mental front in the rapidly growing field of Brownian-time and iterated processes 
and their connections to both (1) new and important higher order deterministic 
and stochastic PDEs and (2) some equivalent fractional Cauchy problems (e.g., 
Allouba and Allouba et al. IH [3 [TJ [2j EJ HH1 M and followup papers; as well 
as the recent articles by Nane (22] [23] , Baeumer et al. [12] , Meerschaert et al. [21] , 
Orsingher et al. [Ml US]; and DeBlassie [S]). I n this new Brownian-time Brown- 
ian sheet PDEs direction, the BTBS-PDEs connections (|1.5[) and (|1.4[) introduce 
a way for a random field to shed light on nonlinear PDEs — and their correspond- 
ing system of n linear PDEs — that are fourth order and interacting; and these 
PDEs, in turn, could be used to gain insight into the underlying Brownian-time 
random field. Guided by our one-time-parameter (n — 1) work in [S] [21 [31 G3 Q] , 
we explore ramifications of such a program in planned upcoming work. In addi- 
tion, following our work with Zheng in |5j, where we initiated the link between 
Brownian-time processes and fractional time derivatives by way of the BTP half- 
derivative generator — an implicit connection to temporally-fractional PDEs that 
was later made explicit and extended in several directions in the works of Nane, 
Meerschaert, Baeumer, Vellaisamy, Orsingher, and Beghin [531 (HI HH (Ml H3] — we 
intend to connect our BTBS with systems of fractional interacting PDEs and their 
corresponding nonlinear fractional coupled PDE in future work. 

1.2.1. Three Brownian sheet-PDEs links. Here, we give the connections between 
the n-parameter M. d - valued Brownian sheet W x and (1) a system of n-linear second 
order PDEs, (2) a corresponding nonlinear second order PDE, and (3) a linear 2n-th 
order PDE. This Brownian sheet-PDEs connection, which characterizes Brownian 
sheet similar to the way the heat equation characterizes Brownian motion, is one 
of the needed ingredients in our proof of Theorem 11.11 and it contains within it 
the first element causing the coupling between u and 'W 3 in Theorem 11.11 above. 
Although some of its aspects are perhaps known as "folklore" among the specialists 
(we obtained the Brownian sheet-PDEs connection — as part of our unpublished 
preprint [7] — and we circulated it as early as 1996), this connection does not seem 
to be recorded explicitly anywhere in the standard literature. We present it here, 
along with a proof that we give in the Appendix for completeness. From a math- 
ematical finance point of view, the system (j 1 . T[) below is at the heart of our new 
Black-Scholes-Merton system for semi-SPDEs bonds and stocks models that we 
introduce and study in a separate joint paper. Due to the fundamental nature of 
Brownian sheet and its connection to space-time white noise, the Brownian sheet- 
PDE connection in Lemma li~T1 below (as well as Corollarv ll . H and Corollarv ll.2p — a 
connection which is as of yet really untapped — and natural generalization thereof 
are useful to further study the interactions between random fields and SPDEs so- 
lutions on one hand and nonlinear PDEs or higher order linear PDEs and their 
corresponding system of linear PDEs on the other (more on that in [7]). 
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The following result gives the Brownian sheet system of PDEs. 

Lemma 1.1 (Brownian sheet linear second order PDEs system). Let {W x (t);t € 
K™} be an n-parameter M. d -valued Brownian sheet on (fi, J^,P) with f{W x (t) — 
x} = 1, for t € <9R™. Let u(t,x) = E,p[f(W x (t))], for a continuous bounded real- 
valued function f : Mr — ► R. Then, u(t,x) is a bounded solution to the system of 
heat-type PDEs 



\\ II u ) AxU > ieN n ,te W'lr- 



;1 V) , (o) |^ = I tt I I '-1 A >" : .i--''h,.t : F:'i ..,■ = 



(&) u(t,x) = /(x); te<9R r |,x€ 

Moreover, u(t,x) is the unique bounded solution to 



(1.8) 



(9-j; / 1 ___ \ O 



9t. I 2 

V i6N„\{j} 



[(6)u v (0,x) = /(x); veR^- x ,xe 



/or each j € N n . In particular, u(t,x) is t/ie unique bounded solution to the system 

(EZD. 



As immediate corollaries to Lemma 11.11 we connect the Brownian sheet W x to 
both a nonlinear second order PDE and a higher order linear PDE. 

Corollary 1.1 (The Brownian sheet second order nonlinear PDE). Under the 
same conditions in Lemma \l.l\ u(t, x) = M,p[f(W x (t))] is a bounded solution to the 
nonlinear second order PDE 

,, 9 , /(S3-(£)-<^<^ "fc."". 

[u(t,x) = /(x); tedWl,xeR d , 

Corollary 1.2 (The Brownian sheet 2n-th order linear PDE). Under the same 
conditions in Lemma \l.l\ u(t, x) = Ep[f(W x (t))] is a bounded solution to the 2n-th 
order linear PDE 

{QTL 

u (t,x) = /(x); t£dW±,xem d , 

where, for each t, the 2n-th order operator S£ n {t, ■), acting on u, is obtained by using 
the Brownian sheet PDE system (|1 .7[) and iteratively applying the tj derivatives for 
j € N„ . In particular, if n = 1 , . . . , 4, then 

JSfi(t,u) = -Au 



(1.11) 



^{t^^^Au + ^A^ 

J%(t, «) = -A« + -t 1 t 2 t 3 A 2 u + -^-^-A 3 u 
z 4 o 

18 5 t 3 t 3 t 3 

~^4(t, u) = -Au + -t 1 t 2 t 3 A 2 u + -tl4tlA 3 u + -±^-A 4 u 
2 4 8 16 
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Remark 1.1. (a) The uniqueness assertion in Lemma 11.11 is not needed in 

our proof of Theorem 11.11 but we give a uniqueness proof in Appendix 
IA1 that judiciously exploits the relationship between Brownian sheets and 
their corresponding space-time white noises to essentially reduce the proof 
to stochastic analysis in one time parameter. 

(b) The uniqueness for (|1.8|) for each j G N n together with the fact that u(t, x) 
solves (|1.7p means that, in the Brownian-sheet case, the system (|1.7| is 
overdetermined or reducible. I.e., reducing the system (|1.7[) to (|1.8[) (for 
any j) leaves the bounded-solutions set unchanged. 

(c) The Brownian sheet heat-type system (jTTTJ) is symmetric in t = (t\, . . . , t n ): 
if 7r(t) is any one of the n\ permutations of t = (ti, . . . ,t n ), then u(t,x) 
and u n (t,x) = u(ir(t),x) satisfy (jl.8|l for each j e N n , without change. 
This is not surprising since the unique solution u(t,x) = Ep[f(W x (t))] is 
symmetric in t. Restricting our attention to t-symmetric bounded solutions 
s(t, x) to the nonlinear PDE (11.91) . we obtain that s(t, x) must satisfy (11.8[) 
for each j, and hence s is the unique solution s(t,x) = Kp[f(W x (t))]. A 
similar observation holds for the higher order PDE (|1.10|) . We leave the 
details to the interested reader. 

Unlike the Brownian motion, "running" the Brownian sheet W x solves heat 
PDEs with variable conductivity. The system of heat-type PDEs (|1.7[) is one in 
which the j'-th PDE (where time is given by tj) may be thought of as modeling 
heat flow in a d-dimensional medium whose thermal diffusivity (conductivity) is 
determined by the product of all other n — 1 parameters ti, i ^ j, in the driving 
Brownian sheet. I.e., the conductivity of the d-dimensional medium varies with 
n — 1 outside factors represented by the parameters £,-, i ^ j, in the driving Brow- 
nian sheet. Probabilistically, this is a manifestation of the fact that i9jVar(W t £C ) = 
IlieN \{j\ ti — which interacts with the Brownian times \B X \ , . . . , \B n \ in the BTBS 
WL, „ (ll) (t) to result in the coupling between u and tyi in (|1.4p and (|1.5|) . Nat- 
urally, when n = 1, the higher order PDE (|1.10[) . the nonlinear PDE (|1.9j) . and the 
corresponding system (|1.7j) reduce to the standard heat equation associated with 
Brownian motion. 

1.3. A new intimately connected Kuramoto-Sivashinsky- variant interact- 
ing PDEs system. Generalizing our KS links to variants of BTPs in |3, 6, 2[[T], we 
modify the BTBS setting above to obtain solutions to interesting systems of novel 
interacting variants of the celebrated KS equation. For a sample of some deter- 
ministic KS results, background, and its applications to turbulence phenomena in 
chemistry and combustion the reader may consult [5S] and the references therein. 
In what follows i = \/~- 1 and should not be confused with the index variable i. 
In our second main result (Theorem II. 2[) ; we show that an interaction similar to 
that obtained for the BTBS-PDEs is also possible in the Kuramoto-Sivashinsky 
setting, we give the corresponding system of linear Kuramoto-Sivashinsky-variant 
interacting PDEs in (|1.16[) . and we give an explicit representations for the three- 
components-solution to each PDE in terms of the complex multiparameter KS 
kernel 



« / n \ n 

(1.12) <« := / exp i £ sAp^ (x,y) "[[K^ds 
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which we call the "transition density" (or kernel) of the n-parameter d-dimensional 
linear Kuramoto-Sivashinsky sheet (LKSS). The kernel 

(1.13) P$(x,y) = - V — ^; x,y£R d , s = ( Sl , . . . ,a n ) £ R n . 

[27nn <= i*i] 
is the n-parameter generalization of the d-dimensional famous propagator (complex 
Gaussian density on R d ) in [3J, associated with Schrodinger equation, and Kf^ s . is 

the density of the BM B l for each iGN„. The kernel K t . x naturally generalizes 

the important linear Kuramoto-Sivashinsky kernel in [3j [6j [2l [1] (the case n = 1 
for d > 1) associated with the linear KS PDE. Working with the complex variant 
of the BTBS density in p. 121) . we now briefly provide a link to such deterministic 
KS-variant interacting system. 

Motivated by the BTBS case above and our KS link in [3,, we define 

v(s, x) := exp I i V" sA / f{y)p^' (x, y)dy 

„ n 

u(t,x)= / v(s,x)l[K™ Si ds 

Jr " \ieN n \{j} J i=i 

® U) (t,x)= [ I J] S A v(s,x)f[K™ St ds 



(1.14) 



for (t,x) e 



£ d and j £ N„. To define u, P, and ^^ on 91 



: IK" we 

adopt the following notational convenience: for t £ R™ and for 7 7^ and 7 c N„ 
we denote by £/(t) the point in W}_ obtained from t = (t\, . . . , £„) by setting £& = 
for every k £ I. For any function g on R™ x R d , g^(t,x) := £/(£/(t),x). With this 
in mind, we give the following temporal boundary definitions: 

(1.15) 



u C/ (t,x) = /(x)exp 



~ E **) 



/CN„ 



V$(t,x) = 



Wg\t,x) 




n w * 

fc6N„\{j} 



exp 




/ = {.?}, 

/ C N n \ {j} , 



n (**-^) 

fe£N„\{j} 



exp 



V km n \{j} J 



.0; IcN n \{j}, 

for I ^ 4>. Of course, if 7 = N„ in (|1.15|l then 7 C = (f> and m^ 7 (t, x) = / (x). Also, if 
n = 1 then ^ c ( / } (t, a;) = U^(t, a;) = u Cl (t, x) = /(x). 
We are now ready to state our second main result. 
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Theorem 1.2 (Linear KS-variant interacting PDEs system). Let f £ C^(R d ;R) 
such that the partial derivatives dijf are Holder continuous with exponent < a < 
1, for alll<i,j <d. Ifu(t,x), U (j) (t,a;) ; and <%W)(t,x) are given by (fLMj) then 
the family 

solves the system of linear Kuramoto-Sivashinsky-variant interacting PDEs 
(1.16) J-u = -^A 2 ^« - ^AU« - \u; t e R^ , x e R d , j e N nj 

O'Cj o Z Z 

with temporal boundary conditions given by (|1.15|) . 

If ri = 1, then ^W) — uu) = U) and (jl.l6[) becomes the linear Kuramoto- 
Sivashinsky PDE (equation (3) in [3]). 

2. Proof of Theorem 11.11 



2.1. A technical lemma: the BTBS case. We start with a diffcrcntiating- 
under-the-integral type lemma that is needed in the proof of Theorem 11.11 



Lemma 2.1. Let W x be an n-parameter d-dimensional Brownian sheet starting at 
x € R d under P; and let f : M. d — > R be bounded and measurable such that d?,f is 
Holder continuous, with exponent < a < 1, for 1 < k, I < d. Fix j £ N„ and let 

(2.1) S^\t,x):=f I J] S , 2 E P [/(r( S ))]n^ s ,rf S 

o 

/or j £ N n , (t, x) € R™ x R d ; where s = (si, . . . , s n ) and Kf^ s . is the transition 
density ofB^(U). Then A 2 <^(t, x) is finite and 

(2.2) A 2 x ^\t,x)= f ( J] S M A^E P [/(^(s))]n^ B ^ds 

,/R + VeN„\{i} / l=1 

If we additionally assume that d^f are bounded (for all 1 < k,l < d), then 
A£<f Ci) (t, a;) is continuous on R'| x R d . 



Proof. For notational simplicity we show that 
Q±g{I) 



(2.3) 



dx 



u) r I \ f) 4 n 

— = / n s u ttt^ if(w x (m n«sjo,a>; * = i,. . .,<*, 

fe jR + \*6N„\{j} / OTfe i=l 

the mixed derivatives cases follow the same steps. In the remainder of the proof, C 
denotes a constant that depends only on n, d, and a and that may change its value 
from line to line. Fix an arbitrary k € {1, . . . d}, and let 

-\x-y\ 2 ' 



exp . 

(2 - 4) ff£&"° ■■= ir,nw*(s) e dv] = - ^i 11 ^;; 



dl/ ^ ' (ZTrnLr^)" 72 
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be the transition density for the Brownian sheet W x , starting on <9R" from x £ M. d , 

o 

under P, and moving to y € M. d at s = (si, . . . , s„) € R" . Using the boundedness 
on / and Problem 3.1 p. 254 in [50] (the case R" x R d with n,d > 1 is a simple 

extension when / is bounded), the symmetry of K s - X )y' in x and y, and the facts 
that 

l im f( v )J^ K BSM = lim — f(v)— K BSin ' d) =0 
(since / is bounded and -grrf{y) is Lipschitz in yk), we get 



d< 



E r f(W*( S ))=^J(y)-^K B ^dy) 



Rewriting the last term in (|2.5|l . and letting /ifc(y) := d 2 f(y)/dy 2 ., we have 

\ -<V 2 



2 -n 



(»* - ;/fc) 2 - FEU »A jrrff 



n 9 



niii 



hk(y)dy 



(2.6) 



El 
Ep 



(^-^( S )) 2 -nr^^ 

n™ s ? 
(^-^(s)) 2 -nr = i^ 

niii*? 



MW x (s)) 

(fc*F(s))-A t («)) 



where we used the fact that Ep Uxk — W^ k (s)) 2 — Yl7=i s = to obtain the last 
equality. Now, using the Brownian sheet scaling, we have 



E T 



(2.7) 



(x k -W^(s)f-l[ S , 



J[s 2 E f 



2 



Wg(s) 

vnr=i s * 



II«?Ef «(!,.• .,!))-! =C[] 



for some constant C. Using (|2.6[) and Cauchy-Schwarz inequality yield 
^E,/(W»(s)) 

/ 2 \ V 2 

(2-8) <(E P f^^^)l!^nkil) Ep 1^(^(8)) -m*)i 2 ] 



< 



mil*? 

-"(&\w'(B)-z\**y'=— 



c 



1/2 



c 



l-a/2 ' 
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where the next to last inequality follows from (|2.7|) and the Holder condition on ht- 
Fixing an arbitrary j £ {1, . . . n}, we therefore have 




A_E P [/(W*(s))] 



(2.9) 



< 



C 



v 2; -v i UL s 



1+3 



II K & t / ^^ds 3 dsy 



«eN„\{j} 



en 



i£N„\{j} t i 



< oo; 



t, 



t e K+ , < a < 1. 



So, (12.3[) . as well as the continuity assertion for d^S^ jdx\, follow by a standard 
classical argument. □ 



2.2. Proof of the BTBS-PDEs system link. Using (JT7J in Lemma O (the 
proof of Lemma ll.ll is given in the Appendix [A"t . Lemma |2. 11 and generalizing our 
BTP-PDEs approach in [5J |3] , we are now in a position to prove Theorem 11.11 



Proof of Theorem \1.1\ It is enough to prove Theorem 11.11 (i) since a solution 
of (|1.4[) is automatically a solution of (|1.5[) . First, fix an arbitrary j e {1, . . . , n} 
and observe that owing to the independence of W x and {i?W} a s well as the 
independence of all the B^'s and the boundedness of / together with Fubini's 
theorem we have 



(2.10) u(t,x) 



f 



' BW,...,B^ 



(t))j=2"/ E r lf(W*(sj)}Y[K™ s Js 



fy {3 \t,x) = 



(2.11) 




n \b {v, {u)\\ f(w^ 1) _ B 

i£N„\{j} / 



w(t) 



n Si E P [/(^)]n^™, Si * 

+ V«eN„\{j} / *=i 



for (t,z) e R^ x R d , where if™ Si is the transition density of B^(t). Differenti- 
ating (|2.10[) with respect to tj and putting the derivative under the integral, which 
is easily justified by the dominated convergence theorem, then using the fact that 
Kf^ s . satisfies the heat equation 

(2.12) 



Of 



BM 









K BM 
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wc have 



(2.13) 



dU 



l (t,x) = 2 n [ E P Jf(W s )]^-K?™ Sj J] K™ Bt da 



*eN„\{j} 



2- 1 / E Ps [f(W s )} ^K™ 8j J] K™ 8i ds 



ieN n \{j} 

Integrating by parts twice in Sj and observing that the boundary terms always 
vanish at oo (as Sj /* oo) and that {d/dsj)Kf^ s .| s =o = but If™ > 0, and 
then using the Brownian-sheet PDEs connection in Lemma 11.11 (II. 7[) successively 
twice we obtain 

J- M (t,x) = 2"- 1 / |^E P [/(V^]nif™ Si ds 

+ 2- 1 / /fBM q ( 5^ [f{w:)] \ jj ***,,*, 

= 2"- 1 / (i J] sA AlE P [f(W*)]f[K?™ Si d S 

jR + \ i€N„\{j} / i=l 

I a\ n si\^m n k ™sA? 

jR + ' \ tefi n \{j} J ieN„\{j} 

n 



in- 1 



V^ 



2 n-3 A 2 



= A 



n * 



v^ 



^ n -A n ^ n ^x^ y 

E + VieN„\{j} / ieN„\{i} 



1 



n ^ 

tg 4 W ' A{3} A g /(x) + lA^^(t,x), 



where pulling A^. outside the integral in the next to last step is justified by Lemma [2TT1 
under the Holder continuity and boundedness assumptions on the second partial 
derivatives of /. The boundary conditions (b), (c), and (d) follow easily from our 
assumptions on / and from the independence of the fiW's and from their indepen- 
dence of the Brownian sheet W x . □ 



3. Proof of Theorem 11.2 



3.1. A technical lemma: the Kuramoto-Sivashinsky- variant case. As in 

the proof of Theorem ll.l[ we need the following differentiating- under-the- integral 
type lemma for the proof of Theorem 11.21 
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Lemma 3.1. Assume f : M. d — > R is bounded and measurable such that d^f is 
Holder continuous, with exponent < a < 1, /or 1 < k,l < d. Let v, 'W^ 1 ' , and 
U^ 6e rfe/?nerf as in (fL~TI)) : 



^«(t,£) := 

(3.1) 

U (j) (t,ir) := 



/ f n ^W^n^u* 

f I 1] Bi)v(B,x)f[Kf$ tBt da 



for j G N„, (t,jc) G M!^ x K d ; and s = (si,...,S n ) G K"; and w/iere i^?£ s . 
zs irae transition density of a one dimensional BM, B^'(ti), starting at 0. Then 
A^(j)(t,x) and A K U^)(t,x) are finite and 



A^«(t,x):= 
(3.2) 

A x U«(t,x):= 



/ f n ^»(«.*)nc* 

jRn \ten n \{ 3 } J i=i 



// we additionally assume that d^f are bounded (for all 1 < k,l < d), then 
A^W(t,x) and A K U (j) (t,x) are continuous on W± x M d . 

In light of our adaptation of the proof of Lemma 2.1 in [3] to that of Lemma 2.1 
in [3] in the one parameter case (n = 1), the proof of Lemma l3.1l is a straightforward 
and similar adaptation of Lemma |2. II above with only minor changes, and we omit 
these details. 



3.2. Proof of the Kuramoto-Sivashinsky-variant interacting system link. 

We now have the needed ingredients for the proof of Theorem II. 2 [ which we present 
next. 



Proof of Theorem [7H Let u, v, U^\ and °^ {j ) be as given in (jl.141) . Differ- 
entiating u(t,a:) with respect to t and putting the derivative under the integral, 
which is easily justified by the dominated convergence theorem, then using (|2.12[) 
and integrating by parts twice using the fact that the boundary terms vanish at 
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±00, we obtain 



J-n(M)=/ ,,(■,*)£**£ [] K™ Si ds 

Ul .i JR n al .i ,,„, v r,i 



*eN„\{j} 






*eN„\{j} 






n*. 



BM 



iGN„ 



(3.3) 



m n *) a >m 

VeN n \{j} / 



Y[ Si ] A x w(s,x) -v(s,x) 
«eN n \{j} 
2 



rfs 



" *eN„ \*eN„\{i} / 

-~A X / II • K '™ -« ( II sAv(s,x)ds~ f ft K™ Si v(s,x)ds 

jR " iGN„ \'»eN„\{j} / " iGN„ 

where we have used the fact that, for each j e N n , v satisfies 



dv 

ds~ 2 



- Yl Si \ A x v(s,x)+iv(s,x) 
\ien n \{j} J 
(3-4) 2 

-4 = -- Y[ S A &lv(s,x)-\ Y[ Si\ A x v(s,x)-v(s,x), 
Sj \im n \{j} J V eN >A{j} / 

and the conditions on / and Lemma 13.11 to take the applications of the derivatives 
outside the integrals in (|3.3[) and (|3.4[) . The temporal boundary conditions are 
trivially satisfied in light of (j!.15[) . The proof is complete. □ 
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Appendix A. Proof of Lemma [TTT] 



Existence proof. First, let K t . x ' be the transition density of an n-parameter 
Revalued Brownian sheet W x (t) = W x (t 1 ...,t„), starting on <9R r | from x E R d , 

o 

under P, and moving to y E R d at t = (ti, . . . , t n ) E R" : 

K ns{n,d) ._ i PIr , i _ i CXD f -\ x -y\' 
K *™ - d y r[W * e dyl - ( 2 7rti • • • t n yn exp { 2 t 1 - t n 



Then, it is easy to see that for any 1 < j < n, 

_d_ K BS(n.d) = 1 

dti t;x ' v 2 



(A] _d_ K BS(n,d)_l A BS(n,d) TT 

iGN„\{j} 

Now, by the definition of u, the assumptions on /, the dominated convergence 
theorem, and (|A.1|) . we readily get: 

II ti ) r f(y)A x K^dy 
^eN„\M / J -°° 

n *i]A, r mK^' d) dy =i( n **W«; 

v i€N„\{j} / - 7 - 00 \J6N„\{j} / 

o 

for j = 1, • • • , n, t E R" . By definition and by the continuity and boundedness of 
/ we easily have lim s ->t u(s, y) = w(t, x) = f(x), for t E 9R" and x E R d , proving 

the existence part of Lemma 11.11 D 



Uniqueness proof. Fix an arbitrary j E {1, . . . , n} and fix arbitrary points t = 
(*i, . . . , tj, . . . , tn) e R™ and x E R d . Let A be Lebesgue measure on <% (R™ _1 ) , let 
@= { J 4cR^ _1 ;A€^'(R™ _1 ),A( J 4) < oo}. Now, let 

(A.2) W j ' x = h^' x (A) = (wjf 1 (A),..., W t j '* d (A)) ; < t < oo, A E 9t\ 

be the d-dimensional space-time white noise associated with an n-parameter Re- 
valued Brownian-sheet W x on a probability space (fi, J?,P): 

Vi? (®*eN»\{,-}[0,ti]) = W&sj); 0< Sj < oo, 

(see Walsh [25] for the case d = 1) where ® denotes the cartesian product of sets, 
and where the j-th parameter Sj is treated as time and all the other parameters 
t? = (ti)ien n \{j} G R™^ 1 are called space (not to be confused with the spatial 
variable x E R d in the PDEs (TO)) , which is the starting point in R d of W x ). Let 
\ &¥ > be the filtration generated by W^ x \ i.e., 

jr« = a ({#;f (A) ; o< 8j <t,Ae #}) ; t > o. 
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As is well known, the space-time white noise W^ x is a continuous orthogonal 
martingale measure |26j . In particular, for any arbitrary fixed ty £ WfT and 
x — (xi, . . . , Xd) £ M d , we have that 



W*M = Wi? (<8« 6 N n \{i}[°.*.-]) '^i ° ^ ^ 



< oo 



is a d-dimensional martingale, in Sj, with coordinates quadratic covariations given 

by 



(A.3) 



w x k %A-),w?n- 



SkiSj Jl U; Sj 6R+, k,l £ {1,. 

iSN„\{j} 



,d}, 



where Ski is the Kronecker delta (6m =1 if k = l and <5;y = if k ^ I). 

Now, fix arbitrary points ty£ M™ _1 and x £ R d . Our first step is to prove 

,i.s M tl s i) '■= u t f (tj ~ Sj,Wt(Sj)), J^ is local martingale in Sj on [0,t,), 

(A. 4) . . 

for any fixed but arbitrary tj > 0, whenever u solves part (a) in (jl.81) . 

To see (|A.4[) . we suppose that u solves the j-th PDE in the system (|1.7j) (a) ((a) in 
(|1.8jl ): and we apply Ito's formula to the function ut y , in the Sj variable, to get 



(A.5) 



ut/tj - Sj , WUsj)) - iH,(tj,W?J0)) = / -djutjt, - r^Wlir^drj 



Vu^-r^WUr^-dWZir,) 



n * 

ieN n \{j} 



AutJtj -r^WUr^drj; a.s. 



where dj is the partial derivative with respect to the j-th variable, where we used 
the dot product notation for the second stochastic integral, and where the third 
term on the right of (JA.5I) is seen by (IA.3I) . This, of course, means that M* (sj) :— 

Ut y (tj — 8j,Wt (sj)), -jp s ] is a local martingale in Sj on [0,tj), since the first and 
the third terms on the right of (|A.5|) vanish because u is assumed to solve the j-th 
PDE in the system lfj~7)) (a). So, (|A.4[) is proved. 
The last step is to show that 

(A.6) If u is bounded and satisfies flES), then u(t, x) = E P [f(W x {t))). 

If u is bounded and satisfies (|1.8j) (a); then, using (|A.4|) . MtJsj), < $j < tj, is a 



bounded martingale. This combined with the assumption that u satisfies {T 
give us that the limit as Sj /* tj exists and must be 

M^(t 3 ) := lim M- y { S j) = f(W?/tj)), a.s. 

s j / *3 

Since M% Isj) is uniformly integrable then 



(b) 



M*( Sj )=E ¥ \MUtj)\^ 



E F 



f{w?M)\*W 



a.s. < Sj < tj. 



Taking Sj = 0, taking expectation, and remembering the notational conventions in 
Notation II .11 we get that 



(A.7) 



(t, x) = u t Jtj,x) = E P Mf (0) = E P f(W*{tj)) = E P [f(W x (t))} , 
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and the proof is complete since j € {1, . . . , n}, t — (ii, . . . , tj, . . . , t n ) G K+, and 
x G IR d are arbitrary. D 



Appendix B. Frequent acronyms and notations key 
I. Acronyms 

(1) BM: Brownian motion. 

(2) BTBM: Brownian-timc Brownian motion. 

(3) BTBS: Brownian-time Brownian sheet. 

(4) BTP: Brownian-time process. 

(5) KS: Kuramoto-Sivashinsky. 

(6) LKSS: linear Kuramoto-Sivashinsky sheet. 

II. Notations 

(1) N n = {l,...,n}. 

(2) Kf^Q s .: The density of a one-dimensional BM starting at 0. 

(3) K t . x ™' : The density (or kernel) of an n-parameter d-dimensional 
Brownian sheet. 

(4) K t . x : The complex kernel or density of an n-paramete <i-dimensional 
linear Kuramoto-Sivashinsky sheet. 
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